Theorem on subwavelength imaging with arrays of discrete sources 
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A theorem on subwavelength imaging with arrays of discrete sources is formulated. This theorem 
is analogous to the Kotelnikov (also named Nyquist-Shannon) sampling theorem as it represents 
the field at an arbitrary point of space in terms of the same field taken at discrete points and 
. . . imposes similar limitations on the accuracy of the image. A physical realization of an imaging 

0^ ' system operating exactly on the resolution limit enforced by the theorem is outlined. 

o ■ 

O ■ PACS numbers: 41.20.-q, 42.25.-p 

(N 

O |. In last decade, certain interest has appeared to electromagnetic systems that, under specific conditions, are able 
to produce electromagnetic fields localized in areas of characteristic dimensions much less than the wavelength of the 
used radiation. Different physical mechanisms can be applied to achieve this, ranging from slabs of Veselago media 
and systems analogous to them P, 0, H, Q to impedance grids or arrays 0, Q or even nonlinear systems 0, [§] . 

The mentioned systems deal with quickly decaying near fields of a source, i.e., with the evanescent modes of the 
• source field. Merlin et al. proposed to transform a part of the propagating spectrum to the evanescent spectrum 

with the help of a planar membrane that modulates the incident field amplitude and phase. The membrane can be 
^ I designed in such a way that the evanescent waves produced by it form a "beam" of sub-wavelength dimensions at a 
O , certain distance from the array. Merlin named this phenomenon radiationless interference. Physically, interference 
c/5 ' here means superimposing the fields of the secondary sources associated with the illuminated membrane. In Q a 
. ^ ■ special distribution of these sources was considered with which the subwavelength field localization can be achieved. 
■ Another interference-based method of subwavelength focusing was proposed in 11, 13, In this approach the 
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incident field penetrates trough a number of slots in a metal screen. The beam patterns produced by displaced slot 
elements form a set of basis functions in which the image-plane field distribution can be expanded and the necessary 
beam formation can be achieved. 

In this letter we will study limitations on subwavelength imaging imposed by discrete nature of the considered 
sources. We will establish a theorem in the space domain that can be seen as an analogy of the well-known Kotelnikov 
(Nyquist-Shannon) sampling theorem in the time domain 3, We will also outline a possible realization of an 



fS| , imaging system operating exactly on the resolution limit enforced by this theorem. 

' Let us start with considering an infinite planar array of high-frequency line currents in free space. We introduce 

' a Cartesian coordinate system with the x-axis oriented along the currents and the z-axis perpendicular to the array 

■rj" ' plane. In these coordinates the currents are at points z = 0, y = nl, where I is the period of the array and n £ Z. 

] The complex amplitudes of the currents are yet unknown, let us denote them /„. Then the total electric field at the 

' point (0,2/, z) produced by this array of currents reads: 
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where k = 27r/A, rj — ^ /io/^o; the time dependence is of form exp(+jwt). 

Let us impose the following condition on the electric field at the plane z — L (image plane): 

Ex{ml,L) = 0, 771 7^ 0. ^ ' 

The respective array and the image plane are shown in Fig. [TJ Physically, Eqs. ^ mean that we would like to 
concentrate the field at the image plane around point y = and at the same time have it completely vanishing at 
discrete points y = ml, m ^ 0. Next, we will look for such a distribution of currents that satisfies ([2l). 
Considering a Fourier series with the discrete currents /„ as its coefficients one can write 

In = ^ f J{q)e-=^Uq, (3) 
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FIG. 1: (color online) An array of line currents at z=0 with the image plane at z = L. We demand the electric field Ex, 
vanishing at discrete points z = L, y = ml, m 7^ 0. 

where J^{q) is yet unknown function. Substituting ([3]) into ([1]) we obtain 

, +00 I 

E,{ml,L) = --^ y H^^\k^{n-mfP + L-^) J{q)e-^''^dq 
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Next, using the Poisson summation formula 
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J2 H^^'iky/n'^P + L2)e-J"« dq = 2j ^ 
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/C(<?), 



(5) 



where /C((7) is a quickly converging series; the branch of the square root in ([5]) is such that V— T = +j. It is seen that 
IC{—q) — IC{q) and IC{q + 27rn) — IC{q). 
From ([4]) and ([5|) we get 



7r 

E,,{ml,L) = -^ j ^(g)/C(g)e-J™« dg. 



(6) 



From ^ and the condition ([2]) it immediately follows that 

Finally, sustituting ([7]) into ([3]) we get the complex amplitudes of the currents: 



4i?o ^ , . ^ , , 1 f cos nq 
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where we have introduced the function 5'(n) that we will call elementary wavelet. In terms of it, condition ^ can be 
rewritten as 



+ C50 



^ '^/{n)H^o^\k^/{n - m)^P + L^) = 

n— — oo 

Now we can formulate a couple of lemmas. 



1, rn = 0, 
0, m ^ 0. 



(9) 
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Lemma 1 For any given electric field distribution E^iyn) defined at discrete points yn = nl in the image plane there 
exists a distribution of line currents defined at the points yn = nl in the source plane that produces exactly the given 
electric field distribution. This distribution of currents is 

4 ^ 

In = -Y H E,{ml)^{n-m). (10) 

m— — oo 



This lemma is obvious: it directly follows from the superposition principle and the solution of the problem considered 
above. 

Lemma 2 — m) form a complete system of linearly independent functions on an infinite set of discrete points so 
that any given distribution of line currents /„ can be uniquely expanded into wavelet series 

oo 

In^ C,n^in-m). (11) 

m— — oo 



The proof of this lemma is given in Appendix. 

Now from the above lemmas we can conclude that the distribution of the field produced by a planar periodic array 
of discrete sources is completely determined by the values of the same field taken at corresponding discrete points in 
any plane parallel to the array plane. Indeed, knowing the discrete values of the field one can reconstruct the currents 
in the source plane with the help of Lemma [1] and then calculate the field at any point of space using ^ . Lemma [2] 
makes sure that the source distribution found in this way is unique. 

Having said that we arrive at the following theorem which is in a sense analogous to the Kotelnikov sampling 
theorem. 

Theorem 1 The electric field of the considered planar periodic array of line currents is completely determined by the 
values of the same field taken at certain discrete points in any plane parallel to the array plane and can be represented 
by the following formula: 

+00 +00 

E,{y,z)= E,i^l,L) Y -^{n-m)Hl;^\k^{nl-yY + z^). (12) 

m— — oo n— ~oo 

As said above, this formula is the result of substitution of Eq. pO)) into Eq. ([T]). 

From this theorem we see that it is impossible to produce arbitrary near-field distributions with an array of currents 
of a given period. The set of realizable distributions is determined by \12\ and is, basically, a set of linear combinations 
of the fields created by the wavelets of currents (fTT|) . Thus, the structure period sets a limit on the amount of details 
that one can expect to observe in a near-field image. Although this fact is rather physically intuitive and has been 
mentioned in the literature before, the above theorem makes it explicit in mathematical terms. 

Let us now discuss some properties of the introduced wavelet functions. The elementary wavelet is given by ([8]). 
When the period of the array is small compared to the wavelength and the parameter L/l is large enough, the series 
}C{q) converges very rapidly, so that only a few first terms of the series should be taken into account. For instance, 
for kl < tt/2 and L/l > 3 it is enough to take only 3 terms of the series. 

Example plots of the normalized elementary wavelets for two cases with different values of the parameters are given 
in Fig.[2l[a, b). One can see that the normalized ^(n) is practically purely imaginary, as the real part of it approaches 
zero for kl < tt/2 and L/l > 3. On the plot there are also shown the distribution of current magnitudes of the source 
introduced by Merlin [9|: cx 1/(1 + n^P/L^) and the distribution that we can obtain from our formulas when 
only a single term of the series ^ with to = is taken into account. We conclude that when L/l 1 the source 
proposed by Merlin can be seen as a one-term approximation of our formulas. This approximation is valid only for 
small n because \E'(n) decays quicker (exponentially) when n — > oo. 

The field distribution in the image plane is depicted in Fig. [2l^c). This function reminds of sinc(y//), however for 
large arguments it decays quicker (exponentially). One can see that the width of the near-field peak at zero level is 
21, as expected. 

Next we would like to consider how the introduced wavelets (which are certain distributions of line currents) can 
be physically realized. It is natural to represent a two-dimensional array of line currents with a grid of thin wires 
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FIG. 2: (color online) (a, b) The normalized elementary wavelets ^{n) as functions of n. The parameters are: a) L/l = 10, 
kl = O.Itt; b) L/l = 5, kl — 0.457r. (c) The normalized electric field in the image plane as a function of y/l. The parameters 
are as in Fig. [IJb). (d) The second term of the total wire impedance normalized to fc77/4 as a function of m. 



illuminated by an incident plane wave. As we would like to excite currents with alternating phases on the neighboring 
wires (see Fig. [Ha,b)) it is natural to load the wires with alternating reactive loads, so that the effective wire 
impedances and wire currents will change from wire to wire in the necessary way. Contrary to the scheme proposed 
in we would like to separate the creation of the phase profile of the currents from the amplitude profile. Indeed, 
the characteristic scale of the amplitude profile L can be much larger than the phase change scale I. In the example 
given in Fig. [S^a) these scales differ by 10 times, so that the phase scale can be subwavelength while the amplitude 
scale is not. Therefore, it is more practical to realize the amplitude profile with the conventional optical means by 
shaping the incident beam with diaphragms and/or lenses, while the phase profile can be realized with wire loading. 

As a proof of concept let us calculate the interaction field and the local field in a wire array with currents distributed 
proportionally to ^'(ti): /„ = /o[5'(7i)/^(0)] for the case depicted in Fig.[2l^b). The interaction field at the m-th wire 
is 

Ei-\ml,0) = -^ ^{nX'\\n-m\kl). (13) 

Given the incident field E}^'^{y,z) the local field acting on the m-th wire can be expressed as E]^^ = E™^{ml,0) + 
E^^^{ml, 0). As we discussed above, we realize the necessary amplitude profile of the currents by shaping the incident 
field accordingly, so that E™'^{ml,0) = E™'^\'^{m)\/'i'{0). Therefore, the total effective impedance (including the 
self-impedance and loading) per unit length of the m-th wire has to be 

_E^_ Eirin^ _ _Jr^ ^ ^MH^^Uln - mlkl) (14) 
~ Im ~ lo *M 4vl/(m) ^^"^^0 ""1^'^- ^^^> 

The phase of the first term in ()14p is alternating from wire to wire while its magnitude remains constant. Moreover, 
if Iq is chosen il7] to be in phase with E™'^ then this impedance term is practically purely reactive (imaginary). The 
dependence of the second term on wire index is plotted in Fig. [DJd) . One can see that the real part of this term is 
constant and equals the radiation resistance of a single wire. Indeed, for the self-impedance of a thin wire of radius r 
we have ^ iZ^^u/iki]) = H^^\kr) w 1 + j(2/7r) log[2/(7fcr)], where log(7) w 0.5772. 

Another important property of it seen from Fig. EJb) is that the imaginary part is also practically constant 
with m [isj . Moreover, for the chosen parameters it happens to be very close to zero. This means that in a real array 
of wires we can neglect the small difference in this part of the impedance and load all wires by loads of only two kinds, 
to account for alternating reactance of the constant magnitude in the first term of . When compared to [lO] this 
greatly simplifies the structure. Moreover, our loaded wire grid is periodic with the period 21, therefore if the beam 
illuminating the array moves along y-axis by an integral number of periods, the subwavelength spot produced by the 
array will also move by the same number of periods. This important property is missing in the realizations proposed 

in [ia,[ii|,[i3,[i3. 
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To conclude, in this letter we have formulated a theorem on subwavelength imaging which is in a sense analogous 
to the Kotelnikov sampling theorem. In the process of derivation of this theorem we solved a problem about a 
distribution of currents in a two-dimensional array that produces a given subwavelegth near-field spot. The found 
distribution of line currents represents an elementary source wavelet. It has been shown that the source proposed in 
[gl. [loj can be understood as a one-term approximation of the elementary wavelet derived in this letter. We have also 
shown that if the necessary amplitude profile of the currents is realized by shaping the incident field by conventional 
optical means the remaining phase profile can be produced with a grid of thin wires loaded by reactive loads of only 
two kinds. When compared to structures considered by other authors such a periodically loaded wire grid appears to 
be much simpler to realize. 



Appendix 

We start with the proof of linear independence. Let us assume that there exists a non-trivial linear combination of 
wavelets such that 

oo 

a^^{n-m) = 0, \fneZ. (15) 

m — — oo 

This linear combination defines a source distribution with all currents equal to zero: /„ = 0, therefore, the field exited 
by these vanishing currents must also vanish everywhere. 

From the other hand, (jl5[) is a non-trivial linear combination, therefore there exists ap ^ 0. But then from 
the definition of ^'(n) it immediately follows that the electric field in the image plane at the point j/p = yl is 
ExiUp) = —{ki]/4)ap 7^ 0. We have arrived at a contradiction, therefore '^{n — m) are linearly independent. 

Completeness of the system of functions — m) can be seen from the following. Consider the linear combination 
of wavelets ([Tl]) with coefficients given by Crn = H^o^\ky/ {p - m)H^ + L^). Thcu, denoting s — n — m and using ([9|) 
we get for the currents 

In= H^^\N{s-{n-pWP + L^)^{s) = \'^^ ^ - (16) 

s— — oo ^ 

Since p is arbitrary we see that the current localized at any point yp = pi can be expressed in terms of the wavelets, 
therefore any distribution of such currents can be expressed also. 
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